Exact Solution of the Isovector Proton Neutron Pairing Hamiltonian 
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The complete exact solution of the T = 1 neutron-proton pairing Hamiltonian is presented in 
the context of the SO(5) Richardson-Gaudin model with non-degenerate single-particle levels and 
including isospin-symmetry breaking terms. The power of the method is illustrated with a numerical 
calculation for 64 Ge for ap/ + gg/2 model space which is out of reach of modern shell-model codes. 
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Exactly solvable models (ESM) provide important in- 
sights into the structure of many-body quantum systems. 
The two main advantages of ESMs are: (1) They can de- 
scribe in an analytical or exact numerical way a wide 
variety of elementary phenomena. (2) They can be and 
have been used as a testing ground for various many-body 
approaches. 

A particular class of ESMs, extensively used in nu- 
clear physics, are the dynamical-symmetry models. In 
this case the Hamiltonian can be expressed in terms of 
Casimir operators of a chain of nested algebras. An ex- 
ample often used to introduce nuclear superconductivity 
(see e.g. Ref. Q) is the rank-1 (Lie) algebra SU(2). Ex- 
amples of dynamical-symmetry models associated with a 
rank- 2 algebra are Elliott's SU(3) model of nuclear defor- 
mation 2] and the SO (5) model of T = 1 isovector pairing 
between neutrons and protons [j| which has found many 
applications in nuclei (see e.g. Ref. 3). 

The concept of quantum integrability, closely linked 
with exact solvability, goes beyond the limits of the 
dynamical-symmetry approach. A quantum system is 
integrable if there exist as many commuting Hermitian 
operators (integrals of motion) as quantum degrees of 
freedom [jj. The set of Casimir operators of a chain of 
nested algebras satisfies this condition. 

Dynamical-symmetry models are usually defined for 
degenerate single-particle levels. Lifting this degeneracy 
breaks the dynamical symmetry but may still preserve 
integrability. The pairing model with non-degenerate 
single-particle levels, of which an exact solution was 
found by Richardson in the sixties @, represents an ex- 
ample of an ESM with such characteristics. Recently, 
more general exactly solvable pairing models, both for 
fermions and for bosons, called Richardson-Gaudin (RG) 
models, have been proposed 0,H|- 

The RG pairing models are based on rank-1 algebras: 
SU(2) for fermions and SU(1,1) for bosons. In this Letter 
we carry out the first step in extending the RG models 
to higher-rank algebras by considering a RG model based 
on the rank-2 algebra SO (5). The model Hamiltonian de- 
scribes a two-component system consisting of neutrons 
and protons interacting through an isovector (T = 1) 
pairing force and distributed over non-degenerate orbits. 



This neutron-proton (np) pairing Hamiltonian with non- 
degenerate orbits has been studied by Richardson |j who 
proposed an exact solution. However, it was shown sub- 
sequently that Richardson's solution is incorrect for more 
than two nucleon pairs [T(i| by explicitly solving the case 
of three- nucleon pairs. Independently, Links et al. de- 
rived an exact solution for the isospin invariant SO (5) 
model by making use of the quantum inverse scattering 
method [ill ■ 

We present here the most general exact solution of the 
RG SO (5) model including isospin symmetry breaking 
terms and for states with arbitray seniority. In addi- 
tion to the construction of the complete set of integrals 
of motion from which more general exactly solvable pn- 
pairing Hamiltonians can be derived, we present here the 
first numerical exact solution of the SO (5) RG model for 
64 Ge in a Hilbert space built from the pf + g g / 2 shells, of 
which the dimension goes well beyond the limits of mod- 
ern shell-model codes based on exact diagonalization. 

SO (5) has also been proposed as the symmetry under- 
lying high-T c superconductivity [l3|. The exactly solv- 
able RG model discussed in this Letter may conceiv- 
ablybe used to generalize SO (5) condensed-matter mod- 
els |l4j by the explicit addition of non-degenerate single- 
particle symmetry-breaking terms. Other possible appli- 
cations might be found in polarized ultracold Fermi gases 
with p-wave pairing interactions [15J . 

We begin by introducing the 10 generators of the 
SO (5) algebra in a representation well suited for nu- 
clear physics problems. Let us define first the three 



T = 1 pair-creation operators 



{h\p\ + p\h\)/\/2, and b\ li = p\p\, where n and p re- 
fer to neutrons and protons, respectively, and i labels a 
single-particle basis (with % its time-reversed state) which 
may be associated with the spherical shell-model basis 
i = jm or with an axially-symmctric deformed basis 
i = am. The three pair-annihilation operators are 6— 1,», 
bo i) an d &+i %■ The three components of the isospin op- 



erator [f +ti = (plhi+pih T )/\/2, f . 
(h\hi + n|rW)/2, and T_ ;i = {h\p t + n\ 
SUt(2) subalgebra of SO (5). These 9 operators together 
with the number operator Ni = p\pi + p\p% + n\fii + n|?% 



i\pi)/y/2] close the 
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define the SO (5) algebra. 

For a system with L single-particle states i = 1, . 
there are L integrals of motion: 



L 1 

Ri = (2 + A)Hi + AT 0) i + 2g £ — (1) 

Zl Zl ' 



where iZj = iVj/2 — 1. The expression Q follows from 
the integrals of motion valid for any semi-simple alge- 
bra of arbitrary rank Since SO(5) is of rank 2, its 
Cartan subalgebra contains two elements namely Hi and 
Toj in the chosen basis which therefore appear linearly 
in QJ. The set of L parameters Zi together with the 
two constants g and A can be freely chosen and it is 
straightforward to check that the integrability condition 
[Ri, Rj] = is valid for any choice of the L + 2 param- 
eters. A simplified version of Q was previously derived 
using the algebraic Bethe ansatz 12]. 

The eigenvalues of the integrals of motion are 
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where e a and tap are solutions of the equations 
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The meaning of the quantum numbers appearing in 
and is as follows: is the seniority of each i level, 
i.e. the number of fermions not paired in time-reversed 
states with isospin T — 1, ti is the isospin of the unpaired 
fermions (this quantum number is often called reduced 
isospin j t = J2i M is the number of T = 1 time- 
reversed pairs, and To is the z component of the total 
isospin, i.e. the eigenvalue of the operator To = J2i ^o,i- 



The total number of nucleons is N = N p + N n = 2M + 
J2i v i whereas their difference is N p — N n = 2T . The 
quantum numbers M, To, Vi, and ti are conserved; T is 
also conserved if A = 0. 

Although any function of the Ri can be used as an in- 
tegrable Hamiltonian, the linear combination X)i=i z i^i 
yields simple expressions for the np-pairing Hamiltonian 



and its corresponding eigenvalues: 
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where C is a constant operator depending on the con- 
served quantities. We have introduced the variables 
£j = Zj/2 and specialized to a spherical basis i = jm. 
The second term on the r.h.s. of breaks the isospin 
symmetry. For A / the operator T 2 does not com- 
mute with the Hamiltonian @ and, consequently, T is 
not a good quantum number but T is still a conserved 
quantity. The same linear combination of the r% and the 
use of the Richardson equations © , yield the eigenvalues 
of©: 
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M+T +t 
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E *"/J + I> [y( 2 + A)-A*j 
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+ |T (T -1). 



(5) 



Each solution of the equations 10 gives an eigenstate 
of the np-pairing Hamiltonian. The spectral parameters 
e a are interpreted as pair energies as in the case of SU(2) 
pairing. However, due to the larger rank of SO (5), a 
new set of spectral parameters wp appears in the equa- 
tions (J3J. These new parameters are associated with the 
SUt(2) isospin subalgebra and there are M + Tq + t of 
them. In the limit A = the number of finite wp pa- 
rameters reduces to M + t — T for each possible isospin 
T. The Bethe ansatz for the SO(5) eigenstates of the RG 
model is a product wave function |l8j : 



M \ M+T +t / M f \ 



0=1 



wp - e a 



where the spectral dependence of the operators is 



(0) 



and J +;Q is a raising operator for bt(e a ): 6j,(e a )J +j0 , = 
^+i( e «) for A* = -1)0, &+i(e Q )/ + , Q = 0, and |A) is a 
lowest-weight state defined by & M ,i|A) = TL^A) = 0. To 
show the behavior of the spectral parameters e a and up 
as a function of the isospin-breaking term A, we plot in 
Fig. n some selected solutions of the Richardson equa- 
tions Q for a system of two neutrons and two protons in 
two shells (jo = 1/2 and j\ =3/2) within the seniority-0 
subspace (M = 2, T = 0, ^ =v t = 0). 

For A = there are two finite cop complex conjugate 
parameters for T = while the two pair energies are real 
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FIG. 1: Pair (dashed lines) and spectral parame- 

ters uip (solid lines) as a function of the parameter A for three 
different states of a 2n-2p system with two shells (jo = 1/2 
and ji — 3/2) in the seniority-0 subspace. The label T refers 
to the isospin in the limit A = 0. The figure shows the lowest- 
energy states of the Hamiltonian @ for each T with energies 
plotted in the bottom-right panel. Since the spectral param- 
eters of the T — 1 state are real, their imaginary part is not 
shown. The interaction strength and single-particle energies 
are g — —1, eo = 0, and ei = 1. 
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FIG. 2: Complex-plane representation of the pair energies e a 
and spectral parameters u)p for the lowest-energy states with 
isospin T = 0, 1, 2 in 64 Ge. The left panel corresponds to weak 
coupling g — —0.05 and the right panel to strong coupling 
g = —0.5. The squares represent the three lowest single- 
particle energies (2e/ = 0.00, 2e P3/2 = 12, 2e fs/2 = 12.5), 
the black circles are the pair energies e a , and the grey circles 
are the parameters up. All energies are in MeV. 



and negative, for T = 1 there is one real, finite u>p and 
two real pair energies e a , and, finally, the T — 2 case re- 
duces to SU(2) for like particles with two e Q parameters 
forming a complex conjugate pair. For finite A there is 
isospin mixing and the number of finite parameters top is 
always two. Figure^thus confirms that the number of fi- 
nite up spectral parameters reduces from M = 2 to M—T 
when A — > 0. For T = 1 one of the real top goes to — oo in 
this limit, vanishing from the Richardson equations PJl 
but giving a finite contribution to the Hamiltonian eigen- 
values (@J 16]. Analogously, in the T — 2 case the two 
u> parameters tend to oo in the A = limit. Also shown 
are the energies of the three eigenstates of the np-pairing 
Hamiltonian (J3J. We emphasize that, while these are the 
eigenvalues of a particular Hamiltonian, the spectral pa- 
rameters completely define the eigenfunction © of the L 
integrals of motion from which H is constructed and 
their corresponding eigenvalues . 

We now turn to the discussion of a numerical calcu- 
lation for 64 Ge. We consider a model space that is well 
beyond modern shell-model capabilities based on exact 
diagonalization: 12 valence neutrons and 12 valence pro- 
tons with a 40 Ca core. The adopted single-particle ener- 
gies are (in MeV) e/ 7/2 = 0.00, e P:i/2 = 6.00, e/ 5/2 = 6.25, 
£ pi/2 = 7.1, and s gg/2 = 9.60, and two pairing strengths, 
g = —0.05 (weak) and —0.5 (strong), are considered. 
We assume isospin symmetry (A = 0) and consider the 



seniority-0 subspace. 

Results for the lowest T = 0, 1, and 2 states are shown 
in Fig. |21 The T = solution corresponds to the ground 
state while the T = 1 and T = 2 solutions are excited 
states in 64 Ge. As in SU(2), the different configurations 
can be classified in the weak-coupling limit. At weak 
coupling (g — > 0) eight pairs occupy the f 7 / 2 level and 
four pairs are in the P3/2 level for the state with T = 0. 
This is reflected in the upper left panel of figure [5] where 
8 pair energies appear close to 2e f T and 4 pair energies 
are close to 2e P3/2 making the corresponding terms in Q 
dominant. Due to the Pauli principle, this configuration 
is not allowed for a state with T = 1 and, correspond- 
ingly, one pair energy is close to 2e/ 5 , 2 . In all cases the 
w/3 parameters are intertwined with the pair energies e a . 
The number of wp parameters (M — T), together with the 
initial configuration at weak coupling, defines each eigen- 
state of the np-pairing Hamiltonian. As \g\ increases, 
the e a and wp parameters expand in the complex plain. 
The solutions are subject to numerical instabilities due 
to singularities arising when a real pair energy e a crosses 
a single-particle energy or when real e a and wp parame- 
ters cross. An example of the first class of crossings can 
be observed in Fig. [3 for T = 2 where the pair energy 
above the p 3 / 2 level at weak coupling goes down with 
increasing g and crosses the P3/2 single-particle energy. 
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FIG. 3: Eigenenergies and occupation probablities of single- 
particle levels of the T = 0, 1, and 2 states in 64 Ge as a 
function of the pairing strength g. 

The T = 1 case shows an exchange of positions on the 
real axis of a pair energy e a and a wp parameter as an ex- 
ample of the second class of singularities. The first class 
of singularities was already present in SU(2) pairing and 
precluded the practical use of Richardson's solution for 
a long time. Recently, a new method to overcome this 
numerical problem was proposed [Tsj . We believe that 
the same procedure can be used to treat the second class 
of singularities as well, allowing the exact solution of the 
SO (5) model for very large systems. 



As a further illustration of the method we show in 
Fig. the eigenenergies and the occupation probabli- 
ties of single-particle levels as a function the pairing 
strength. The occupation probabilities can be obtained 
making use of Hellmann-Feynman theorem which ex- 
presses them in terms of derivatives of the eigenvalues 
of the integrals of motion r*j as (N Pt i) = 1 + and 
(N nti ) = n - - (1 + A) (1 + |g) + 1. These deriva- 
tives can be related to the derivatives of the spectral pa- 
rameters e a and oj a , which in turn can be obtained taking 
the derivatives of the Richardson equations 10 ■ 

In summary, as an application of generalized RG mod- 
els, we have presented the complete solution of the SO (5) 
isovector np-pairing problem. The generalization allows 
the introduction of one-body symmetry-breaking terms, 
such as non-degenerate single-particle energies, yielding 
an exact solution of the SO (5) np-pairing model for ar- 
bitrary seniorities even if it includes an isospin-breaking 
term. The numerical solution of the SO (5) Richardson 
equations was presented for the specific example of 64 Ge, 
together with a discussion of the behavior of the spec- 
tral parameters for weak and strong pairing. With this 
work the exact solution for large systems with SO (5) sym- 
metry is now available which could be of great impor- 
tance in condensed-matter physics when addressing the 
phenomenon of high-T c superconductivity 0, 0] . Fi- 
nally, the treatment of higher-rank algebras like Sp(6) 
and SO (8) opens the possibility of exact nuclear struc- 
ture calculations with more realistic quantum integrable 
models. 
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